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ABSTRACT 



We study the radial ionization structure at the surface of an X-ray illuminated 
accretion disk. We plot the expected iron Ka line energy as a function of the Ed- 
dington ratio and of the distance of the emitting matter from the central source, for 
a non-rotating and a maximally-rotating black hole. We compare the predicted disk 
line energies with those measured in an archival sample of active galactic nuclei ob- 
served with Chandra, XMM-Newton and Suzaku, and discuss whether the line energies 
are consistent with the radial distances inferred from reverberation studies. We also 
suggest using rapidly- variable iron Ka lines to estimate the viscosity parameter of an 
accretion disk. There is a forbidden region in the line energy versus Eddington ratio 
plane, at low Eddington ratios, where an accretion disk cannot produce highly-ionized 
iron Ka lines. If such emission is observed in low-Eddington-ratio sources, it is either 
coming from a highly-ionized outflow, or is a blue-shifted component from fast-moving 
neutral matter. 
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1 INTRODUCTION 

Recent high-resolution spectral studies by Chandra, XMM- 
Newton and Suzaku have revealed that emission lines at 
6.7 or 6.9 keV (probably emitted by highly ionized iron) 
are common features in the X-ray spectra o f active galac- 
tic nuclei (AGNs, e.g., 



Pounds et ail 120031: iReeves et all 
20041: lYaqoob fc Padmanabhanl 12004 iBianchi et all 1200a 



Shu et al-lhoiCtlPatrick et al.ll20ldl )7The Ka emission lines 



produced by ionized iron in opticall y thin, photoionized ma - 
terial in AGNs have been studied bv lBianchi fc MattJ (2002). 
Their results may require an iron overabundance by a factor 
of a few to account for the observed line strength. 

The energy and ionization balance at the surface of an 
accretion disk can change due to the X-ray illumination, 
and this affects the line emission. Thus, it is important to 
determine the radial ionization structure of X-ray photoion- 
ized accretion disks, in order to study the reflected spectra 
and the associated iron Ka emission. The energy of the iron 
Ka line generally increases as iron becomes m ore stripped, 
from 6.4 keV for Fe I to 6.9 7 keV for Fe xxvi l|Kaspi et al.1 
l2002l : iPaerels fc Kahnll2003l ). This provides a diagnostic of 
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the accretion disk structure: we can compare the calculated 
photoionization structure of the disk with the observed en- 
ergy of the iron Ka line, assuming that the line was emitted 
from the disk. 

In this letter, we calculate the radial ionization struc- 
ture of an X-ray illuminated AGN accretion disk as a func- 
tion of Eddington ratio, AEdd = L^ i/ L^dd, where Lboi is 
the bolometric luminosity and I/Edd is the Eddington lumi- 
nosity. We then compare the predicted range of ionization 
parameters and iron Ka line energies with a sample of re- 
cent line observations. We plot the predicted and observed 
line energies E a as a function of AEdd- This gives us clues on 
the origin of the iron Ka lines. Throughout this Letter, we 
assume standard cosmological parameters ho = 70 km s _1 
Mpc~\ n m = 0.27, and ft A = 0.73. 



2 METHODS AND RESULTS 

X-ray reflection from the surface of cold matter 
around compact ac creting obj e cts h a s been studied 
by m a ny authors dBasko et alj 1 19741: Guilbert fc Reesl 



19881: iLightman fc White! Il988l; IWhite et~ 



George fc Fabianlll99ll ; IMatt et al1ll99ll; iDone et al.lll992l: 
Ross fc Fabianl Il993t iMatt et al.1 1 19931 : ICzernv fc Zvckil 



1988; 
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19941; iKrolik et all Il994l; iMagdziarz fc Zdziarskil Il995l: 



Ross et al 
J— 



Reynolds fc Begelman 19971; Blackmar] 
20001: lYoung fc Reynold j 



iBallantvne et al 



2000 



2002 



1996; . 
19991: iNavakshin et all 
Navakshin fc KallmanT 120011 

Miniutti fc Fabian 20041 ; etc). Following those studies, we 
assume an idealized lamppost geometry for the X-ray illumi- 
nation of the reprocessing matter: all the X-rays are emitted 
isotropically by a point source lo cated on the symmetry axi s 
above a standard accretion disc (|Shakura fc Sunvaevll 197*3 ) . 
at a certain height h. Ignoring the effects of light bending in 
the vicinity of the black hole (BH), the ionization parameter 
is 



£ = 8.97 x 



f(r)r 



-3/2 



(1) 



T)* J - ■ ( r 2 +^2)3/2 

(|Matt et al.lll993l ; iRoss fc Fabianlll993r i, where r is the ra- 
dial coordinate on the disk, f(r) = 1 — yj&jr, -q is the to- 
tal radiative efficiency, ?y x is the X-ray fraction (L x = rj x 
AEdd-I/Edd), and a is the viscosity parameter. In Eq.(l), h 
and r are in units of r g = GMbh/c 2 . The efficiency n = 0.06 
for a Schwarzschild BH and n = 0.31 f or a maximally spin- 
ning astrophysical BH l|Thorne I I1974T ). Here we assume a 
standard value of viscosity parameter of a ~ 0.1, although 
a large range of values has been used in previous stud ies, 
from a ~ Q.Q1 dMiller fc Stonell2000l ; IStarling et al.ll2004| ) to 



0.3 (|Esin et al.lll997l ) 



We took the values of the iron Ka line energy emit- 
ted by an X-ray i lluminated disk a s a function of £ from 
the calculations of iMatt et al] j 19931 Fig lb). Note that the 
relation between E a and £ is almost independent of the in- 
cident angle of the illuminating flux; here we assume that 
E a is only a function of the iron ionization state. Using 
Eq.(l), we calculated E a as a function of AEdd for differ- 
ent values of h and r (Figure 1). Plotted as black lines are 
the results for a Schwarzschild BH; the blue lines are the 
results for a maximally-spinning Kerr BH. We repeated the 
calculation for six values of r between 5r g and 200r g in the 
Schwarzschild case, and the same six values of r in the Kerr 
case. The radial values are plotted as labels in Figure 1, 
next to the Schwarzschild curves. For each value of r, we 
performed our calculation for three values of h — 5r g , 8r g , 
and 15r g . The curves corresponding to those three values of 
h are plotted as solid, dashed and dot-dashed lines, respec- 
tively. 



3 COMPARISON WITH OBSERVATIONS 

We compared the predicted and measured values of the line 
energy in the (-E Q ,AEdd) plane (Figure 2; only the calcula- 
tions for a Schwarzschild BH are plotted, for simplicity). 

The Eddi ngton ratios A^ dd of the observed sample ar e 

taken from lZhou fc Zhand (|2010h and lZhou fc Wand l|2005h . 
who studied a large sample of X-ray luminous AGNs. Using 
our theoretical curves, we can directly estimate the effective 
distance of the line-emitting region from the illuminating 
X-ray source, in the various sources. For instance, obser- 
vational values falling on the dashed line corresponding to 
(r = 35r g , h = 8r g ) indicate a distance R = (r 2 + /i 2 ) 1//2 « 
36 r g . This distance determines a characteristic time lag 
between primary X-ray source and line emission. From the 
observed values of time lag and intrinsic line width, one can 




^bol/^Edd 

Figure 1. Predicted line energies in the (E a , \^dd ) plane. The 
black lines are the results for a Schwarzschild BH; the blue lines 
are for a maximally spinning astrophysical BH. In each case, we 
repeated the calculation for six values of r (indicated as labels 
next to the Schwarzschild curves), and three values of h = 5r g , 
8r g , and 15r g . For each r value, the three corresponding h curves 
are plotted as solid, dashed and dot-dashed lines, respectively. 




Figure 2. Comparison of the observed Fe Ka line energies with 
the predicted values. For simplicity, here we show only the (r, h) 
families of curves calculated for a Schwarzschild BH. The radial 
distances and the meaning of the solid, dashed and dot-dashed 
lines are the same as in Figure 1. The E a values observed in 
narrow-line Seyfert 1 galaxies (NLS1) are plotted as blue stars; 
those from broad-line Seyfert 1 galaxies (BLS1) are magenta tri- 
angles; those from quasars are orange squares. Multiple line en- 
ergies are plotted for most galaxies. See Table 1 for details of the 
observed values. 



estimate the B H mass and spin, with the iron line r everber- 
ation method (|Revnolds et al.ll 19991 ; iLiu et al.ll201fj|) . In the 
rest of this section, we compare the line emission distances 
inferred from the (-E a ,AEdd) plane with those estimated in 
the published literature, for a sample of Seyfert galaxies. 

NGC7314: simultaneous observations of the NLS1 
NGC7314 with Chandra and RXTE have revealed vari- 
ability in the spect ral features on a timescale < 12.5 ks 
|Yaaoob et al.| [2003). corresponding to a light-crossing dis- 
tance of w 500r g and a Keplerian radius of w 18.5r g for a 
BH mass of 1O 6 ' 7O M l|Zhou fc Wandl2005l ). This is in ex- 
cellent agreement with the Fe xxv emission region (~ 20r g ) 
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estimated for this object from our line energy curves (Figure 
2). 

Mrk 766: the X-ray spectra o f the NLS1 Mrk 766 show a 
broad emission line at « 6.7 keV ([Pounds et alj|2003h . From 
our plots in Figure 2, we infer that the Fe xxv emission 
region is located at 35rv < r < 50r g from the central BH. In 
addition, I Turner et alj |2004l 'l found a transient narrow line 
at « 5.6 keV and interpreted it as evidence for blob ejection 
of neutral or low-ionized material. If so, the rest-frame line 
energy is m 6.4 keV, and the line emission region may be 
located at 50r E < r < 100r g , correspo nding to a distance o f 



10 14 cm for a BH mass of 10 6 6 M (IZhou fe Wan 



in broad agreement with the results of lTurner et al 



2005), 



2004) 



1H0707— 495: the spectrum of this NLS1 galaxy shows 
a sharp drop at e nergies w 7-7.5 keV (|Boller et all [2002: 
iGallo et all l2004bl ). A partial covering model was intro- 
duced to reduce the need for an extreme iron overabundance 
jTanaka et al.ll2~004h . However, an alternative possibility is 
that the X-ray spectrum is dominated by ionized reflec- 
tion rather than ab sorption. In this scenario, it was argued 
(|Fabian et all 120041 ') that the two XM M-Newton observa- 
tions su pport the ligh t bend ing model of iMiniutti fe Fabianl 
( 2001). iFabian et all (|2009h also measured a time lag of 
about 30s between the soft energy band (0.3-1 keV) and 
the medium energy band (1-4 keV). Combining this short 
time lag with the measured width of the broad iron 1 emis- 
sion, they suggested that the X-ray reverberation comes 
from matter very close to the event horizon of a rapidly 
spinning BH. However, consideri ng a larger energy range 
(0.3-7.5 keV). iMiller et all (|2010l ) argued instead that the 
observed time delays extend up to about 1800s in the hard 
band; this is consistent with reverberation caused by scatter- 
ing of X-rays passing through much more dis tant absorbing 
mate rial For our adopted mass of 10 6 - 37 M Q (|Zhou fe Wand 
l2005f ). A Ed d ~ I This suggests that 1H0707-495 may have 
a highly-ionized disk within ~ 50r g for a Schwarzschild BH, 
or within m 20r g for a rapidly spinning BH. In our scenario, 
those characteristic distances may be the origin of the ob- 
served iron Ka lin e s at rest-frame ene rgies « 6.5-6.7 keV 
jFabian et al.ll2009l : IZoghbi et ailbOlOT ). 

NGC3516: simultaneous Chandra an d XMM- 
Newton observations showed (|Turner e t al. 2002) two pairs 
of weak emission features, including a component at ~ 6.9 
keV, symmetrically located around a strong, narrow 6.4 k eV 
line. This structure was interpreted (I Turner et al.ll2002T ) as 
evidence for relativistic broadening of disk lines from three 
different radii. In our model, we find that for AEdd ~ 0.01 
(Figure 2) even the innermost region of the accretion disk 
must remain neutral. Thus, if the observed 6.5 keV and 
6.9 keV lines are emitted from the disk, they must be 
blue-shifted peaks of a 6.4 keV line, rather than being 
emitted by ionized iron. This supports the interpretation of 
iTurner et aH (|2002T ). 

Mrk 841: we infer that the observed 6.4 keV line is emit- 
ted from a region « 50r g from the central Schwarzschild 
BH (or closer, for a Kerr BH). This corresp onds to a time 
lag o f 20 ks for a BH mass of 10 790 M (|Zhou fe Wand 
l2005l l in good agreement with the variability timescale < 36 
ks found in two XMM-Newt on observations l|Petrucci et all 
120021 : lLonginotti et alj[2004l l 



4 DISCUSSION AND CONCLUSIONS 

Iron at the surface of an accretion disk is significantly ion- 
ized when the Eddington ratio AEdd is larger than a critical 
value. Assuming a standard viscosity parameter a ~ 0.1, 
the critical value above which iron in the innermost part of 
the disk becomes ionized is Asdd ~ 0.1 for a Schwarzschild 
BH, and AEdd ~ 0.3 for a maximally-rotating astrophysical 
BH. We studied the radial ionization structure of an X-ray 
illuminated accretion disk, and calculated the energy (in- 
creasing with the ionization parameter) of the iron Ka lines 
emitted from the disk. We plotted those energies as family 
of curves in the (E a , AEdd) plane, parameterized in terms of 
radial distance of the emitters and height of the illuminat- 
ing X-ray source above the disk plane (lamppost model), for 
a non-rotating and a maximally-rotating BH. We compared 
our model with the observed Ka line energies from a sample 
of AGNs. 

A substantial fraction of AGNs show highly-ionized 
iron Ka emission. The origin of the ionized emission is 
still debated. Our results suggest that it may come from 
two different sources: the accretion disk (for AEdd > 0.1) 
or the photoionized material in the outflow (for AEdd 
0.1). Our model presented here is based on simple assump- 
tions, such_as_TOnstant density without vertical stratifica- 
tion (|Matt et al.lll993h . but our main goal is to illustrate 
an important physical effect, which is unlikely to depend 
substantially on the details of the disk structure. 

The critical AEdd depends on a, which parameterizes 
our ignorance of detailed accretion physics (I Ji et all 120061 : 
IMiller et al1l2006l '). Despite forty years of observational, ex- 
perim ental and theoretical studies since IShakura fe Sunvaevl 
<| 19731 ). we are still unable to determine the disk viscosity 
accurately. The theoretical dependence of the observed iron 
Ka line energy on a suggests that we can reverse the argu- 
ment: if we have independent measurements of a BH mass, 
spin and luminosity, we can estimate a using the ionization 
curves in the (-E a ,-^Edd) plane, by combining the informa- 
tion on centroid energy and rapid variability timescale. It is 
plausible that a ~ 0.1 is in agreement of iron line observa- 
tions of a few AGNs. 

Iron near the disk surface cannot be ionized at low 
accretion rates and low Eddington ratios. There is a for- 
bidden region in the (-E Q ,AEdd) plane, below which ion- 
ized Ka line emission cannot come from an irradiated disk. 
Observationally, several low-luminosity AGNs in that re- 
gion show Ka emission features at 6.5-6.9 keV. We argued 
that such features are either coming from a highly-ionized 
outflow, or are blue-shifted components from fast-moving 
neutral matter. Alternatively, the intermediate energy line 
(6.5-6.7 keV) seen in NGC3516, NGC7213 and Fairall 9 
may come from an evaporating/condensating region, as pre- 
dicted by the disk transit i on model at AEdd ~ 0.01-0.02 
ilRozanska fe Czernvl |2000| : iLiu fe Taanl 120091 : iQiao fe Liul 
l2009h . 

Finally, we emphasize that the current observations do 
not yet allow us to put robust constraints on the origin of the 
highly-ionized iron Ka lines. The line parameters derived 
from X-ray spectral fitting are strongly model-dependent, 
and have large uncertainties. Future observations with the 
next generation of X-ray space telescopes (such as the pro- 
posed GRAVITAS mission) will resolve the profile and con- 
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Table 1. Observed rest-frame central energies of the iron Ka lines in our sample of galaxies. Col.(l): common name of the object; 
Col. (2): redshift, from the NASA/IPAC Extragalactic Database (NED); Col. (3): source of the X-ray data; Col. (4-6): centroid energies 
(f fixed); Col. (7): Eddington ratio, from Zhou & Zhang (2010) and Zhou & Wang (2005); Col.(8): references for the X-ray spectra: (1) 
Patrick et al. (2010); (2) Gallo ct al. (2004a); (3) Takahashi ct al. (2010); (4) Schmoll ct al. (2009); (5) Bianchi ct al. (2004); (6) Bianchi 
et al. (2005); (7) Fabian et al. (2009); (8) Turner et al. (2002); (9) Reeves ct al. (2004); (10) Pounds ct al. (2003); (11) Reeves ct al. 
(2001); (12) Reynolds ct al. (2004); (13) Longinotti ct al. (2003); (14) McKcrnan & Yaqoob (2004); (15) Matt ct al. (2001); (16) Pounds 
ct al. (2001); (17) Bianchi et al. (2003); (18) Yaqoob ct al. (2003); (19) Pctrucci ct al. (2002). 
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